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Quantum noise in a model of singly resonant frequency doubling including phase mismatch and 
driving in the harmonic mode is analyzed. The general formulae about the fixed points and their 
stability as well as the squeezing spectra calculated linearizing around such points are given. The 
use of a nonlinear normalization allows to disentangle in the spectra the dynamic response of the 
system from the contributions of the various noisy inputs. A general "reference" model for one- 
mode systems is developed in which the dynamic aspects of the problem are not contaminated by 
static contributions from the noisy inputs. The physical insight gained permits the elaboration of 
general criteria to optimize the noise suppression performance. With respect to the squeezing in the 
fundamental mode the optimum working point is located near the first turning point of the dispersive 
bistability induced by cascading of the second order nonlinear response. The nonlinearities induced 
by conventional crystals appear enough to reach it being the squeezing ultimately limited by the 
escape efficiency of the cavity. In the case of the harmonic mode both, finite phase mismatch and/or 
harmonic mode driving allow for an optimum dynamic response of the system something not possible 
in the standard phase matched Second Harmonic Generation. The squeezing is then limited by the 
losses in the harmonic mode, allowing for very high degrees of squeezing because of the non-resonant 
nature of the mode. This opens the possibility of very high performances using artificial materials 
with resonantly enhanced nonlinearities. It is also shown how it is possible to substantially increase 
the noise reduction and at the same time to more than double the output power for parameters 
corresponding to reported experiments. 

PACS numbers: 42.50.Dv, 42.50.Lc, 42.65.Ky 



I. INTRODUCTION 



Second Harmonic Generation (SHG) has nowadays quite a long tradition as a mean of squeezed light generation 
The preferred experimental setup has been the doubly resonant configuration as, at least in principle, permits 
arbitrarily large squeezing. However, such scheme has been hampered by the technical difficulties arising from keeping 
the resonance in both modes simultaneously. Thus, in spite of the development of very ingenious stabilizing procedures 
||, for the moment it has been only possible to maintain the double resonance for a few seconds. Certainly, this kind 
of experimental delicacy can hardly surprise when dealing with the generation of non-classical states of light. In view 
of such difficulties, some experimental efforts have been recently redirected to singly resonant configurations 
Although, the maximum noise suppression is then limited to a 90% the efforts resulted in very stable intense 
squeezed light sources with degrees of squeezing even surpassing those reported in the doubly resonant counterparts 
|i|. This evolution highlights the importance of reducing to a minimum the technical demanding of new proposals in 
a so experimentally challenging field. 

At the same time, singly resonant Optical Parametric Oscillation (OPO), the most successful method to squeeze the 
vacuum has been generalized to singly resonant Optical Parametric Amplification (OPA), i.e. a laser driving 

in the harmonic mode has been added, again showing an extraordinary stability at quite high noise suppression 
values in the fundamental mode []i"o| . Although the squeezed beams are in this case much less intense than the in 
SHG counterpart, this setup permits a control of the phase of the squeezed quadrature, something which allowed a 
spectacular demonstration using quantum tomography, of the different kinds of squeezed states . 

In view of this experimental success it seems timely to extend the quantum mechanical model beyond the pure 
phase matched cases. More specifically, we address here quantum noise reduction in an extension of the conventional 
singly resonant SHG to include also a coherent input in the harmonic mode as well as phase mismatch between the 
interacting waves. 

On the other hand, an increasingly number of papers is being devoted to study quantum noise in systems combining 
different kinds of nonlinearities (see, for instance, |l^-[l6| for some recent contributions). In particular, the combination 
of x*- 2 "* with Kerr-like nonlinearities in c.w. cavity systems has been quite extensively studied [ p^|Jl7| -^]. Even 
exact full quantum results have been obtained showing, for instance, the emergence of tristability not present in the 
classical counterpart Ejj . With respect to the squeezing performance the results appear as very promising at least 
in degenerate doubly resonant configurations ]15|,p2[ . The simplest system from the implementation point of view, 



1 



combining this two kinds of nonlinearities that the authors can think of is precisely a singly resonant second order 
nonlinear system with phase mismatched interacting waves as then, by virtue of the cascading effect, an effective 
Kerr-like third order nonlinearity appears. 

In order to make the search of strong noise reduction through the parameter space affordable we are bounded to 
the standard linearization procedures, the only capable of yielding analytical results. Inside the linear approximation 
perfect squeezing is possible at dynamic instabilities. We make use of this fact to find optimum working points showing 
up maximum squeezing. They are, however, an artifact of the method as the linear approximation breaks down at 
the instabilities. What matters for the practical implementation of new squeezed light sources (our ultimate goal) 
are the optimum paths through the parameter space approaching such points. Let us explain a little more what we 
mean. Fixed a particular parameter there will be a set of values of the remaining parameters (including the frequency 
as such) tuned up to yield the maximum noise reduction. When the chosen parameter is varied an optimum path is 
defined by the set of parameter values maximizing the noise reduction at any stage of the variation. The real optimum 
working point will be somewhere along these paths before reaching an instability. Thus, these paths would guide the 
experimentalist towards the optimum working point in the real experimental setup. The essence of our approach 
to find such paths will be to isolate the dynamic aspect of the squeezing behavior from the static contributions to 
the noise coming from the different inputs. A simple adequate normalization will disentangle the two aspects of the 
quantum noise behavior. This simplifies the analysis sufficiently to allow a characterization of the optimum paths. 

Another crucial issue when dealing with the squeezing performance of a system is precisely the election of the most 
relevant parameter to compare the different configurations with each other or with the reported related experiments. 
There is no universal criterion to determine the squeezing efficiency of a given device. An efficient setup regarding 
power consumption, i.e., when compared for fixed input power could well be deceptive when compared for the same 
output power and perhaps inadequate to some spectroscopic applications. However, within the state of the art of 
the present squeezed light generators, the main concern is to improve the squeezing figures themselves having other 
considerations such as the power consumption a relative importance. Under this perspective probably the parameter 
of utmost importance as far as c.w. resonant systems are concerned is the energy load inside the cavity. Indeed, the 
usual causes of squeezing degradation such as blue-light-induced red absorption come from an excessive mean photon 
number inside the cavity capable of significantly degrade the material optical response at the relevant frequencies. 
These considerations will lead us to define another normalization this time useful for the evaluation of the squeezing 
efficiency with respect to the intra-cavity photon number. 

The sketch of the article is as follows. In section || the quantum mechanical model is presented. In section [II 
the evolution equations are linearized, the fixed points of the system obtained and their stability studied. Section IV 
gives all the formulae regarding quantum noise spectra in the system. In section [v| a general approach to one-mode 
systems is developed which allows the definition of general criteria to characterize the optimum paths and applied 
to the specific case here addressed. Finally the limits of the model and possible implementations are thoroughly 
discussed in section VI, concluding the article with a summary of the most relevant results obtained. 



II. QUANTUM MECHANICAL MODEL 



The system we want to address consists in a second order nonlinear medium coupling two modes of frequency w 
(fundamental) and 2oj (harmonic) respectively and placed inside a ring cavity resonant only with the fundamental 
mode. We will also assume just one input-output mirror of finite reflectivity. The effect of phase mismatch when only 
the fundamental mode is driven has been experimentally studied in p5| where bistability induced by cascading was 
demonstrated. The classical evolution equation of the fundamental mode, a, as given in p5[ , reads 

^ = -b + i6 + ^(Afc)M 2 ] a + j¥f c a in . (2.1) 

The nonlinear coupling depends on the wave vector mismatch Afc = k(2uj) — k(oj) as K(Ak) = 
2 f Q m /* it* (Afc, z)u(Ak, z')dz'dz/L 2 m being L m the length of the nonlinear medium, u{k,z) the spatial dependence 
of the resonator mode and v is proportional to the second order nonlinear susceptibility (see below). Splitting K in 
its real and imaginary parts, Eq. ( |2.lD can be recast as 

^ = - [7 + [i\a\ 2 + i(8 + T\a\ 2 )] a + ^¥f c a m . (2.2) 

For a plane wave geometry 
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FIG. 1. The dependence of Ki and K T with respect to the phase mismatch. 
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where K r and Ki denote the real and imaginary part of K(Ak). In this way the nonlinear dynamics is divided in a 
nonlinear absorption (the up conversion of photons) and a nonlinear dispersion (the cascading effect). The behavior 
of both parameters with Ak is plotted in Fig. |l|. Notice that for any finite Ak the nonlinear dispersion is also finite, 
periodically completely do minating (null frequency doubling). 

Quantization of Eq. (2^) is then accomplished independently for each effect. Regarding the nonlinear absorption we 
use the two-photon model proposed in [^6| , while nonlinear dispersion is accounted for by a fourth order Hamiltonian, 
H = (hT/2) a^ 2 a 2 as in the standard theory of optical Kerr effect. It represents a Hamiltonian modification of the 
two-photon absorption model so that the quantum mechanical equation reads 



— = — [7 + iS + (/! + ir)a t a] a + 2^fjia)bi n + 



(2.4) 



where Latin characters denote the annihilation operators for the corresponding classical (Greek characters) modes. 
Two extra terms not present in the classical analog appear, namely, a white noise input, Wi n , accounting for the 
fluctuations induced by the scattering and the absorption in the crystal (7,5 — 7 — 7 C ) and a parametric "gain" term 
coming from the, classically empty, incoming harmonic mode, bi„. Eq. (2.4) is complemented with the boundary 
conditions 01 



bout = - b v 



(2.5a) 
(2.5b) 



from which the output spectra can be computed. Input fields are assumed to be in coherent states. In particular, 
allowing a coherent state different from the vacuum for the incoming harmonic mode we generalize the system to the 
case of driving both modes. In the case Y = 5 = 0, the squeezing properties as well as the applicability to quantum 
nondemoliton measurements of this system have been studied in detail in |27j . 

The used definitions for the creation operators give the following relations with the usual experimental param- 
eters (see appendix in Q): the input and output powers are P w ,i„/ „t = hoj{a\ n/out a in / out ) and P 2u) ,in/out = 



h2uj{b\ 



in I out 1 in l 'out I 



the circulating power is %ui{a^a) /r being r the round-trip time and the single-pass power- 



conversion efficiency (in W ) is 2t v/Tiuj. 



III. LINEARIZED EVOLUTION EQUATIONS AND LINEAR STABILITY ANALYSIS 



Defining fluctuation operators as 



a = a + 5a , 

bin,out — Pin, out H~~ ^in,out t 



(3.1a) 
(3.1b) 
(3.1c) 



a linearization of Eqs. (2.4) and (2.5) yields 

[ 7 + i ^ + 2(^ + i r)|a| 2 ] 6a + [2^/Jl[3 m - (// + il> 2 ] fa* 



d5a 
~~dT 



and 



+ 2^fjia*5b in + yj2r/~ c 8a in + V^Ts w ™ > 
<56 ont = 2a^/JISa - 5b in , 



(3.2) 



(3.3a) 
(3.3b) 



being aj niOU t, f3i n . ou t the me an v alues of the corresponding input and output modes and a a stable fixed point of the 
classical counterpart of Eq. (2.4), i.e. 



da 
~dt 



[7 + iS + (p + iT)\a\ 2 ] a + 2^fjia*(3 m + ^/2q c a in . 

Equating to zero the l.h.s. of Eq. ( |3~4| ) a "state equation" for the fixed points is obtained, namely, 

_ \Z2t7 {[7 + fin - i(6 + Tn)] a m + 2^ffjf3 in a* n } 
° ~ ( 7 + /m) 2 + (<5 + rn) 2 -4/i|/3 m | 2 ' 



(3.4) 



(3.5) 



with n = \a\ . Let be 9, <f> and tp the phases of a, a^ n and respectively. Then, dividing both sides of Eq. (3.5) by 



'2 lc 



\ a \ e W-*m [( 7 + ^ n )2 + {5 + Tn f _ 4/ i|/3 m | 2 ] 
(7 + /m - i(£ + Tn)) e^"^ + 2^ |/3 m |e~^~ v/2) 



(3.6) 



Taking the squared modulus in both sides a quintic equation for n is obtained 

= n [(7 + /in) 2 + (5 + Tn) 2 - Afi\l3 ln \ 2 ] 2 - 2 7c |a m | 2 {(7 + /in) 2 + (6 + Tn) 2 + 4/i|/3 m | 2 + 
4VM I An I [(7 + M n ) cos(20 -tp) + (S + Tn) sin(20 - 95)]} . 



(3.7) 



The re al an d the imaginary part of Eq. (3J3) determine the sin(# — tp/2) and cos(# — tp/2) as functions of the solutions 
of Eq. O) 



cos(6> - tp/2) = 
sm(0 - tp/2) = 



\a\ 

\Oin 

\a\ 



'2 7t 

/27c 



(7 + fin + 2y/ji\pi n \) < 



y>/2) + (<5 + Tn) sin(0-y>/2) 



( 7 + /m) 2 + ((5 + rn) 2 -4/i|/3 m | 2 
(7 + an - 2^/71 |/? m |) sin(0 - y/2) - (S + Tn) cos(0 - y>/2) 
(7 + /in) 2 + (<5 + rn) 2 -4/i|/3 m | 2 



(3.8a) 
(3.8b) 



Eq. (3.7) allows for numerical calculation of the fixed points given the input fields. But it can be interpreted also as 
a linear equation for |ct; rl | 2 , i.e. 



27c|ai 



n [(7 + fin) 2 + (S + Tn) 2 - 4/i|/3 m | 2 



| 7 + fin + 2^jl\[3 m \eW-^\ 2 + A^JI \f3 m \(S + Tn) sin(20 - p) ' 



(3.9) 



The positive character of the r.h.s. is not always guaranteed and therefore not for every value of the parameters a real 
positive n is possible. Notice, however, that in the cases in which this happens a simultaneous change of the sign of 5 
and r yields a consistent set of parameter values. As we shall see, this fact will have useful consequences in regarding 
the analysis of the quantum noise behavior in the system. 

The stability of the fixed points is g overned by the real part of the eigenvalues of the drift matrix associated with 
the linearized evolution equation (|3.2[). Very simple algebra yields 
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A± = -(7 + 2/xn) ± ^J\{n + iT)a 2 - 2^M A«| 2 - (<* + 2Tnf . (3.10) 

Provided that the real part of both eigenvalues are negative the fixed point will be stable. With respect to the 
phase-matched SHG case (T = and = 0, always stable), although both T and 8 alone tend to stabilize the 
dynamics, in combination are able of destabilize the system. A finite /3j„, on the other hand, can promote instability 
depending on its relative phase with respect to a, the case 9 — ip/2 — ±ir/2 maximizing the effect. All of these new 
instabilities, however, correspond to zero eigenvalues without a finite imaginary part. In other words, contrary to the 
double resonant SHG there is no Hopf bifurcation and, consequently, no selfpulsing solution. 



IV. SQUEEZING SPECTRA 



For a given quadrature of the electric field, Xg ut (t) = a out (t)e~ 
the noise spectrum of such a quantity, i.e. 



t (i)e l6ls , the squeezing spectrum is simply 



S(u) = C 
= C 



(SXgf {t)5X$f (t + T))e~ u 
(5XS ut (w)5XS ut (-w'))du>' . 



(4.1) 



being C some normalization constant and the averages are assumed stationary. As a function of the annihilation and 
creation operators Eq. (|4.1|) is rewritten as 

(4.2) 



S(u) = C (5al ut (ui)8a out (-uj)) + Re{exp(-i29 s )(6a out (u>)Sa out (-oj))} 



where use has been made of the stationarity of the average and Re denotes real part. From this expression it is evident 
that the noise is minimized and therefore the squeezing effect maximized for a quadrature phase such as 



corresponding to a phase 



S(ui) = C (5al ut (oj)5a ut(-v)) ~ \(5a out (uj)6a out (-u!))\ 



V(bj) — 7T 



(4.3) 



(4.4) 



where v{u>) is the phase of (Sa ou t (m) Sa ou t (—u>)). The spectrum of the conjugate quadrature (i.e. with a phase v{oS)/2) 
corresponds to a plus sign in Eq. (4.3) and by virtue of the Heisenberg princi ple s hows an excess noise above the 
vacuum. Taking C — 1 (corresponding to vacuum noise units) and splitting Eq. (4.3) into a vacuum noise component 
plus a normally ordered part we finally arrive to 



S-,+ (u>) = 1 + (: 5al ut ((jj)8a out (-uj) :) T 10 8a out (uj)Sa out (-uj) :)| 



(4.5) 



for both the squeezing and the "stretching" spectra. After tedious but simple algebra, the spectra of the fundamental 
and second harmonic modes can be written as 



S b _ + (cj) = l + 8/j,n\B 

where B = 2^fp,f3 in — (fi + iT)a 2 and 

7V_,+ = 2|B|(7 + 2^n) qp J + 2fin) 2 - (5 + 2Tn) 2 - 



D ' 
D 



\B\ 2 + oj 2 } 2 + 4( 7 + 2fm) 2 (5 + 2Tn) 2 . 



■ U) 
\2, ,2 



D = [(7 + 2fj,n) 2 + (S + 2Tn) 2 - \B\ 2 - lu 2 ] + 4( 7 + 2\xnfup- 

The correlations defining the squeezing phase v(lo) are given by 

(6a oa t{cj)6a oat (-uj)) = 4 7c B [lo 2 + \B\ 2 + (7 + 2 l in) 2 - (5 + 2Tn) 2 + i2(-/ + 2 f in)(5 + 2Tn)] /D , 

(Sb out (u)8b out (-uj)) = 8/1 a 2 B [uo 2 + \B\ 2 + (7 + 2fin) 2 - (S + 2Tn) 2 +i2(j + 2/m)(<5 + 2Tn)] /D . 

The trigonometric equations for the corresponding phases are quite complicated and rather useless. However, an 
interesting consequence can directly be drawn from Eqs. ( |4.8| ), namely, for detuning-s such as S + 2Tn — the phases 
are independent of ui equaling those of B and a 2 B respectively. 



(4.6a) 
(4.6b) 

(4.7a) 
(4.7b) 

(4.8a) 
(4.8b) 
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V. SQUEEZING PERFORMANCE 



In the previous sections we have developed the general raw formulae regarding quantum noise in the system. This 
section is devoted to the analysis of the quantum noise behavior implied by them an in particular to proceed with our 
program of finding optimum quantum noise reduction. However, before any specific assessment of the quantum noise 
performance we will elaborate a little more on the formulae (mainly by adequate normalizations) in order to gain 
physical insight and ease our task. As an aftermath we shall obtain general results going well beyond the specifics of 
the system addressed here. 



A. General results concerning one-mode systems 



Let us begin defining a nonlinear and a total decay rate as j n i = 2/j,n and 74 = 7 + j n i respectively. We shall 
scale the evolution with this total decay rate defining an dimensionless time r = 74 1. In the spectra (4.6) the only 
dependence on 9 is through B disappearing for (3 in = 0. It is also possible to restrict this dependence to such a term 
directly in Eq. (3.2) and Eq. (3J3) by means of appropriate phase shifts of the modes. All together account for 



ddc 
~dV 



= -\l 



iA] 5c + B 5c< 



(5.1) 



where the tilde represents divided by 74, A = 5 + 2Tn, B = 2^/JlSi,, 



c = 


a e , 






din, out 




Cin,out — 








din, out = 


bin, out 

\/it 


e- l2e 




Win _ 




Sin = 







In agreement with the previous notation 5i n denotes the mean value of 
modes are 

8d out = y/2%i 5c - Sdin ■ 



[fM + iT)n and the modes are redefined as 

(5.2a) 
(5.2b) 

(5.2c) 
(5. 2d) 

The boundary conditions of the new 



(5.3a) 
(5.3b) 

We shall 



For coherent states, the correlations of the new input modes remain as white noise but in the scaled time r 
refer to the previous formul ae a s the tilde normalization. 

The evolution equation (5.1) encodes the dynamic response of the intracavity system to a series of noisy input 
channels (5di n , 5ci n and Sj n ). Quantum Mechanical consistency, i.e., conservation of equal-time conmutators, imposes 
a fluctuation-dissipation relation which under this normalization reads 



7ni + 7c + 7s = 1 ■ 



(5.4) 



The evolution equation (|5.1| ) along with Eqs. ( |5.3| ) are now written in such a way that the input-output couplings are 
real- valued as in the standard input-output formalism [^,^9). This is a completely general result. Provided a well 
defined linearized theory in the sense of preserving equal-time conmutators we only need the adequate set of phase 
shifts of the input channels (a trivial unitary transformation preserving conmutators) making the couplings real- valued 
to obtain a theory formally equal to the standard input-output formalism simply because this is the theory preserving 
the equal-time conmutators when the couplings are real-valued. Thus, for any system with only one effective mode 
there is a formulation in which the intracavity field follows 



dSc 
~dT 



= - \l + iA] 5c + B 5 J 




(5.5) 



with 



G 



JV 

= (5.6) 

71=1 

The frequency scale jt denning the dimensionless time r is just the real part of the factor multiplying 5c after the 
phase shifts. N-l of the input channels will have a time- reversed counterpart corresponding to the outgoing channels 
fulfilling 

8<» lt = ,/vfc5c-8<&. (5.7) 

The remaining input channel will account for the irreversible losses. The corresponding spectra are related with the 
intracavity spectra by 

= 1+ : Sl i+ (u>) : = 1 + 2 7 „ : S_,+ (£) : , (5.8) 

where : S— +(Q) : denotes the intracavity spectra and we have made use of the proportionality of normally ordered 
intracavity and outgoing correlations [ p9| . The spectra S 1 ™ + (<D) coincide with the spectra of the original formulation 
as the new outgoing modes are just a phase sift of the originals. 

Let us now define a sort of "reference" system with only one time-reversible input channel, i.e., 

^ = -[l + iA]5c + B5ct + V25<g f , (5.9) 

and 

5c r l f t =V2Sc-5cZ f . (5.10) 
Obviously : S"x(w) : = 2 : S- t +(u>) : so that we get finally 

S^ + (Q) = l+%:S r _ e f + (Q): . (5.11) 

This is the central result of this section. Let us elaborate a little about its interpretation. Squeezing in a given 
output channel means that for a certain range of phase shifts the corresponding quadratures show an intensity of their 



fluctuations below that of the associated incoming channel (assumed in a coherent state). In view of Eqs. (5.7), the 
amplitude of the outgoing fluctuations is a coherent superposition of the intracavity and the incoming fluctuations. 
Squeezing is possible if an adequate correlation between 5c and the relevant input channel is established. But the 
intracavity field is nothing else than the dynamic response of the intracavity system to the incoming channels. The 
input channels are uncorrelated and so the dynamic response of the intracavity system to them. A given input channel 
can consequently correlate only with the dynamic response to itself. The presence of any other input channel can 
only degrade the effect. The great advantage of the tilde normalization is that makes this fact explicit. Indeed, Eq. 
( |5.11 ) express the output spectra as the dynamic response of the system to an isolated input channel, i.e., : S^ + (Cj) :, 
scaled down by the "static" contribution to the noise owing to the presence of extra input channels. The scale factor 
7„ is just the ratio between the coupling constant of the chosen output channel and the sum of all of them. 



Eq. (5.11) greatly simplifies our task of finding the optimum path to m axim um noise reduction as we can center 
our efforts onto the simple reference system described by Eqs. ( |5.9| ) and ( 5.10) ). Even more interesting the results 



concerning the reference system will be of general applicability to any one-mode system including as such any multiply 
resonant system under adiabatic elimination of all the modes but one. The normally ordered spectra of the reference 
system are easily calculated as 



_ 2|£|t\/(1 + w 2 + |B| 2 - A 2 ) 2 +4A 2 

: S ref ACu) : = A\B\— ^— ^— . (5.12) 

' ' ' (l-w 2 -|B| 2 + A 2 ) 2 +4d. 2 V ' 

Our first step is to determine if the dynamic response is capable of a t otal n oise suppression. Perfect squeezing can 
only occur at a dynamic instability. Equaling to zero the l.h.s. of Eq. ( 3.1 0| ) (the only possible unstable eigenvalue) 



and after proper normalization an equation determining the instability can be written as 

1 + A 2 = |B| 2 . (5-13) 

Written in this way an interesting parallelism with the standard OPO below threshold shows up, i.e. an instability 
appears when the modulus of the "losses" coefficient equals that of the "parametric" coefficient, a sort of natural 
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extension of the condition for the instability in the conventional OPO for which the coefficients are real. Inserting the 



instability condition (5.13) in : S_(u>) : results in 



_ 2\B\ - \Ja\B\ 2 +Cj 2 (u 2 +4) 
w 2 (w 2 +4) 



■ s? f (-) : = m v . _ . (5 . 14) 



Applying L'Hopital's rule with respect to ci> 2 , : Sj £ ^(Q) : equals -1 at ui = 0, that is, perfect squeezing is obtained at 
the instability, again in parallel with OPO. In other words, the dynamic response of the system assuming that the 



condition Eq. (5.13) is reachable, is capable of a complete suppression of quantum noise. 

Spectrum (5.12) is simple enough to permit analytical optimization. Taking partial derivative in : S_(ui) : with 

respect to u> and equaling to zero, ui = appears as the optimum point whatever the values of A and \B\. The same 
applies to A = when taking partial derivative with respect to A. Notice that this last condition implies also a 
squeezing phase independent of the frequency. The optimized noise obtained imposing these two conditions simplifies 
to 

: ^ t: = -(iTW' ( 5 - 15 ) 

with a minimum at the instability |B| = 1 approached monotonically. These conditions (u) = 0, A = and \B\ = 1) 
will help us in finding optimum paths. In particular, moving \B\ from zero to one while maintaining A = uj = 
defines an optimum path reaching the instability for the "reference" model. 

An optimum path is defined solely by the squeezing spectrum leaving aside the "stretching" one. It is important to 
study also the accompanying excess noise on the conjugate quadrature for it could invalidate in practice the optimum 
path if this excess noise is unbearable high. The minimal excess noise production imposed by the Heisenberg principle 
corresponds to S- (lj)S+(lu) = 1. It is a perfect complementary relation between quadratures: the deamplification 
of fluctuations in a given quadrature must equal the amplification of fluctuations in the conjugate. In that case we 
are dealing with a Mini mum Uncertainty State (MUS) for those quadratures, the text-book definition of a squeezed 



state. Adding 1 to Eq. (5.12) and after some minor algebra 



2\B\ ±\J(Co 2 + \B\ 2 + 1 - A 2 ) 2 + 4A 2 

S r _ e {&) = ± = '— . (5.16) 

' + (1-w 2 - |5| 2 + A 2 ) 2 +4£ 2 V ' 

Straightforward algebra leads to S r£i \u>)S , ^ J F (o)) = 1. The excess noise is minimum again in parallel to the standard 
OPO system with real coefficients. 



B. Standard normalization 



As discussed in the introduction we are here principally interested in the squeezing behavior with respect to the 
photon number n. Unfortunately the tilde normalization is inappropriate to such a task as the frequency scale depends 
on n itself. It is far more convenient to use 7 -1 as the time scale instead of 7 t ~ 1 and to normalize the photon number 
as m — i/n/~/. In complete parallelism to the tilde normalization we have then, 



dSc 
"d7 



- 1 + 2mK r 
+ 2 yj mK r 8di n 
where the hat represents divided by 7 and 



i(S + 2mK. 
\/2% 5c 



5c- 



l K r r\i n - (K r + iKi)m 



<5c t 



+ V 2 7s s i 



V* 7 3 e -m 

7 



(5.17) 



(5.18) 



which represents the harmonic mode input amplitude normalized to the value at the standard OPO threshold. The 
spectra (O) become now 



8 



Sl i+ (u>) = 1 + 4%\B\ 



7V_ 



+ 



S°_ :+ (u) = l + 8K r m\B\ 



D 

N-,+ 

b 



with 



N_. + = 2\B\{l + 2K r m) T \j [(1 + 2K r m) 2 - (<S + 2*Qm) 2 + \B\ 2 + w 2 J + 4(1 + 2K r m) 2 (S + 2K l m) 2 , 
D = [(1 + 2if r m) 2 + (<5 + 2K l m) 2 - |S| 2 - lj 2 } 2 + 4(1 + 2K r m) 2 Co 2 , 



(5.19a) 
(5.19b) 

(5.20a) 
(5.20b) 



and B — \J K r rji n — (K r + iKi)m. We shall refer to the above formulae as the hat normalization. To refer the 
frequency to the cavity decay constant with the subsequent re-normalization of the system parameters, i.e. the hat 
normalization, is quite a standard procedure in the literature. 

Some care must be taken when studying the quantum noise behavior as a function of n (or m). It is not a free 
parameter of the problem as would be the input fields or the phase mismatch but it is in a nonlinear relation with 
them. We need, therefore, to check that the proposed values of n are indeed a solution of Eq. fl3.7| ). Fortunately the 
spectra ([O]) do not depend on the overall sign of 6 + 2Tn and therefore the conclusions reached in section III about 
the existence of \cti n \ permit a safe variation of n in search of strong noise reduction provided the stability of the 
corresponding fixed points. 



C. Squeezing at the fundamental mode 



Applying Eq. (5.11) to the fundamental mode 



S* i+ (£>) = l + %:S^.(t5) 



(5.21) 



It is clear that the best performance corresponds to 7„; = 0, that is, either n = or /i = 0, as then j c maximizes 
to rj = 7 C / (7 C + 7 S ) (the escape efficiency of the cavity). The case n = corresponds to the very well known case 
of squeezed vacuum generation. For /i = and finite n the system is formally equivalent to a resonan t opt ical Kerr 
effect system whose quantum noise behavior has been amply studied previously [ fj0| . The condition (|5.13| ) reduces 
for /i = to S = — 2nT ± yj n 2 T 2 — j 2 , the well known turning points of optical dispersive bistability |30| but with 
the nonlinear dis per sion induced by cascading. Indeed, such cascading induced bistability has been experimentally 
demonstrated in |25|. Rewriting it within the hat normalization the condition reads 



6± = -2mKi ± v m 2 K 2 - 1 , 



(5.22) 



where Ki = — \/tt as it is evaluated at AkL m = 2ix. Once 7 C is independent of Ak and n, the optimum path with 
respect to m (the only remaining free parameter) is determined solely by the reference system. It corresponds to 
increase m till m = ir (where the condition (5.22) is reached) while maintaining 5 = 2m/ 'ir and u — 0. Fig. [2] displays 
the evolution of both the maximum squeezing and the maximum excess noise following such a path for three values 
of the escape efficiency, namely, 0.9, 0.99 and the ideal 1. The noise is expressed in dB's with respect to the vacuum 
noise. A Heisenberg limited excess noise appears in such specular image of the squeezing. The instability 

is signaled by the divergence in the excess noise. Above it, the curves shown are not physical as they correspond to 
unstable fixed points. The case 77 = 0.99 in Fig. || shows an excellent behavior with an almost Heisenberg limited 
excess noise till near the instability. 

Fig. U illustrates the idea of optimum path by comparing the 77 = 0.99 plot of Fig. || against various cases with fixed 
values of 6. Below m — n, for a given m the maximum squeezing is obtained when A = as expected. Above m — -k 
it is not possible to reach the minimum noise of 1 — rj fulfilling A = 0. 

Being /1 = the formulae simplify enough for allowing a simple expression for the squee zing phase. More specifically, 
from B = — iTa 2 , 9 m = 6 + tt/2. On the other hand substituting Eq. (3.5) in Eq. ( [2.5a| ) results in \f2rf c a ou t = 
a (7c — 7 S + iTn), giving a squeezing phase relative to that of the output field of 

7r / Tn 

arctan 

2 



7c ~ 7s 
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80 




m 

FIG. 2. Noise spectra at zero frequency of the fundamental mode following and optimum path for three escape efficiencies 
of the cavity including the ideal case 77 = 1. The curves above the divergences are not physical. 
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FIG. 4. Squeezing in the fundamental mode at zero frequency as a function of the phase mismatch for a low intracavity 
photon number. 



At the instability Tn — 7 and for low 7^ it approaches 45°. There is a possibility of an extra control of the squeezing 
phase not present in the conventional Kerr effect system by making use of the harmonic mode. Taking a finite \x 
but low enough so that fin -C 7 and at the same time a /3j n high enough to imply 2^JJi[3i n ~ 7, still we will have 
7„; k, 7 while the squeezing phase relative to the output field will depend on both the modulus and the relative phase 
between the input fields. In practice, however, maintaining fin very low could imply a exceedingly high (3i„ in order 
to have a 2^JJif3i n intense enough for a significant influence on the final phases. Of course, the instability point would 
accordingly depend on 

There is no hope of any behavior similar to the reported in ]l5| as a competition between second and third order 
nonlincaritics needs fi 7^ 0, opening the fundamental mode to the fluctuations of the input harmonic mode with 
strong deleterious effects. At most, some remnants of the enhanced efficiency coming from the competition between 
nonlinearities can be observed for low 7„; . Then, as shown in Fig. ||, the best working point is not necessarily located 
at u = 0, i.e. maximum squeezing is obtained with a finite mismatch. 



D. Squeezing at the harmonic mode 



For the harmonic mode Eq. ( 5.11 ) yields 

S b _ t+ (Q) = l + %i:S r _ e ? + (u): . (5.23) 

Now, the situation is the complete opposite: the performance is favored by a finite fi in order to have a non-zero j n i 
and a large n to approach the ratio j n i to one. In fact, under ideal conditions of perfect dynamic noise suppression 
and no absorption and scattering losses (7^ = 0), the squeezing in both modes are complementary in the sense of 

1 S«+S*=2-^-^ = 1, (5.24) 
7t It 

a direct consequence of the fluctuation-dissipation relation (p.q ). This complementarity has been previously reported 
for the doubly resonant degenerate parametric oscillator |3l|]. The maximum squeezing available for the harmonic 
mode whatever the dynamic response of the system is easily obtained by setting : S^{l.(u>) : to -1 in Eq. (5.23), that 
is, 

2mK r 1 

S M = 1 — = • 5.25) 

This static contribution to the nois e is n ow nonlinear in the sense that it depends on the phase mismatch and m. 
An immediate consequence of Eq. ( [5.25 ) is the possibility of an arbitrarily large quantum noise reduction for any 



finite value of K r . The 1/9 limit of the conventional phase-matched SHG is therefore due to a failure of the setup 
to maximize the dynamic response of the system. Let us center then, firstly in the SHG-like case with /3j„ = as 



11 





0.3 




V.Z 




0.10 






Pi") 





(N 










-0.1 




-0.2 




-0.3 




2k 



4k 
Ak L 



6k 



8k 



FIG. 5. The value of Ki — 3K r as a function of the phase mismatch. 



it includes the above mentioned conventional setup (the experiments in 
given by (directly in the hat normalization) 



and 0). The instability points are now 



-2mKi ± 



yJm2(K?-3K?) 



AK r m - 1 



(5.26) 



Both kinds of nonlinearities (dispersive and absorptive) are in this case necessary as the factor Kf — 3Kf needs to be 
positive to allow S± to be real. The phase-matched case is therefore excluded. 

Fig. p| shows Kf — 3Kf as a function of the phase mismatch and indeed near above n it is positive. Optimum 



approaches to (5.26) are now more difficult to evaluate than in the fundamental mode as both the dynamic processing 
of the noise and the static contribution from the noise inputs (encoded in 7„;) depend on to and Afc. With respect to 
to to is clear that the static part is optimized at to — > oo. This limit can be approached letting Q = and S = —IKim 
(i.e. A = 0). \B\ reduces in this case of (3i n = to m^jKf + Kf /{ X+2mK r ) showing a monotonic increasing behavior 
with respect to m from to the maximum (at m — > oo) 



\B\ = - 
1 1 4 



1+1^ 

' K r 



(5.27) 



Notice that for Kf — 3Kf it consistently equals 1. We have then, both j n i = 1 and the fastest approach to 1 of 
\ B\ w hen m — > oo. The refore , the squeezing along an optimum path with respect to Afc is given by substituting Eq. 
(5.27) in the spectrum ( 5.15] ) and then the obtained : S op t ■ in Eq. ( |5.23| ). Obviously, in a real experiment to can be 
large but always finite. Let us take as a "large" to one giving a Sm around 20 dB as in the case r\ = 0.99 of Fig. ||. 
This corresponds to to = 50. Figure ^ displays 5^ + as a function of the phase mismatch in such a case. To illustrate 

the modulation exerted by Sm we have take this time 5 equal to the real part of Eq. ( j5.26| ) plus a very small number. 
In this way the plot remains valid for the whole range of the AkL rn . While Eq. ( 5.26| ) is complex the condition A = 
is almost fulfilled and above the instability the noise suppression reduction follows Sm- Again the pernicious effect of 
the instability regarding the excess noise has a very short range. For comparison Sm is also depicted. 

The optimum path with respect to m is much more complicated to find because the intricate dependence of Ki and 
K r with respect to the phase mismatch. Figure [7] has been generated finding numerically the minima of 5*^(0) while 
scanning the range of to. For comparison the phase-matched case is also depicted showing an asymptotic behavior 
towards —10 log 9. For low values of to the effect of j n i overwhelms the dynamic response so that the best value 
corresponds to maximize K r . As soon as the two curves depart from each other the dynamic response dominates the 
behavior and the minimum noise is at the instability as in Fig. |^. At this stage the optimum path begins to follow the 
instability all the time. It should be taken then, as a mathematical limit. However, in view of Fig. H, before reaching 
it, bearable values of the excess noise are accessible with a slight diminution of the squeezing. 

It is worth to mention that a squeezing as large as 48% induced by cascading has been very recently reported p2| . 
The cascading was due, however, to a detuning of the pump mode in a triply resonant non-degenerate OPO with a 
much lower finesse for the pump mode rather than by phase mismatch. Under such conditions a cascaded is a lso 
induced leading ideally to perfect squeezing in the pump mode. 
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FIG. 6. Noise spectra at zero frequency of the harmonic mode following a nearly optimum path with respect to the phase 
mismatch for the SHG like case. 
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FIG. 7. Squeezing in the harmonic mode along an optimum path with respect to the normalized intracavity photon number 
(m) for the SHG like case compared with the phase-matched SHG case. 
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FIG. 8. Noise spectra at zero frequency (harmonic mode) following an optimum path with respect to the normalized input 
harmonic amplitude (rji n ). The curves are not physical above the divergences. 



Although a finite mismatch allows to reach Sm > the overall optimum working point corresponding to Ak = is out 



of reach. The question arises then, of if it is possible to fulfill Eq. (5.13) at Afc = 0. A glance at the definition of B 
suggests it should be possible adding a driving to the harmonic mode. In such a case 7„j simplifies to 2m/ (1 + 2m) 
obviously independent of r\i n . Constructing an optimum path with respect to the harmonic input reduces then, to set 
A = Co = 0. Phase matching along with A = implies K r = 1, i<Q = and 5 — so that the instability condition 



(5.13) simplifies to 

1 + 2m = \rnn - m\ , (5.28) 

a perfectly achievable condition. We can further optimize by choosing the phase of rji n adequately to approach B to 
one as much as possible. The extreme cases correspond to r\i n real, i.e., T)i n = — (1 + m) and rji n = 1 + 3to. From Eqs. 



(3.8) and ( |3.9| ) it is easy to check that they correspond respectively to — ip/2 = ir and (f> — ip/2 = 0. The negative 
case maximizes B. It has been previously reported in |27|] . Taking squared modulus of Eq. ( |2.5b|) the negative case 
appears as promoting harmonic output power while the converse is v alid f or the positive. The squeezing phase is also 
easy to calculate in this case. In particular, given the correlation ( 1.8b ), v(lo) is determined by the phase of c?B 
(independent of w as A = 0), something proportional to 

(Vin ~ m)e M . 

The corresponding squeezing phases are 9 m = 29 + it for the negative case while for the positive case it changes from 
9 m = 26 + it to 9 m = 29 + ir/2 at rji n = m. On the other hand, the output harmonic amplitude is proportional to (see 
Eq. (pg )) 



bout oc (r) in - 2m)e'( 28+I ' . (5.29) 

Consequently, the relative squeezing phase for the negative r\i n is 7r, i.e., amplitude squeezing. The positive case 
is more complicated. It remains equal 7r (amplitude squeezing) till 77 = m. Above this value it changes to +ir/2 
depending on the sign of r\i n j2 — m yielding in any case phase squeezing. At a first glance, it appears there is a sudden 
change from amplitude to phase squeezing when the input phases are fix to — (p/2 = and \j]i n \ passes through 
to. It is not so however, as at this point B = and the state collapses to a coherent state with no squeezing. The 
situation is clearly depicted in figure || where 5^. + (0) are displayed as a function of r\i n assumed real. The r.h.s. of 
the plot corresponds to 0— ip/2 = while the l.h.s. to 0— ip/2 = ir and for negative ordinates it should be considered 
as an optimum path with respect to 77^. The behavior is completely symmetric with respect to r\i n = m where both 
the squeezing and the excess noise equal that of the vacuum. 

The optimum path with respect to to is now given by Sm at K r = 1. Figure | is the equivalent to Fig. |t] for the 
new situation. It represents the maximum efficiency as far as quantum noise reduction is concerned the system can 
yield in any way with respect to to. The improvement with respect to the standard phase-matched SHG as well as to 
the optimized SHG is certainly high. 
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FIG. 9. Maximum squeezing (harmonic mode) as a function of the normalized intracavity photon number (m) in nonlinear 
second order singly resonant device. For comparison the phase-matched and optimized SHG cases are also shown. 



VI. DISCUSSION AND CONCLUSIONS 



Two are the main purposes of the present work. On one side, to gain physical insight about the origins of quantum 
noise in singly resonant systems. On the other, to explore their potential as squeezed light sources. In such a task we 
have used a model including all the relevant physics we wanted to address but simple enough to be tractable. The 
results shown in the previous sections certainly reveal a high potential of the studied configurations. An evaluation of 
the limits of the model in reproducing the real physical situation as well as a discussion of possible implementations, 
seems, therefore, in order. 

One obvious idealization of the model is to assume perfectly coherent inputs neglecting the excess noise of real 
lasers something expected to have deleterious effects at low frequencies. White and coworkers |33| have developed 
an analytical approach to this problem resulting in an impressive agreement with the experiments. As expected, the 
excess noise completely destroys the squeezing at low frequencies. In their experiments, however, the deleterious effect 
was restricted to only 7 MHz by adding a mode cleaner to the system, the spectrum coinciding with the ideal one 
out of this range. Even better, in |11| the laser noise was shot-noise-limited down to 1 MHz, again using an external 
mode cleaner. 

Considering as sensible the assumption of coherent states for the input modes as well as a value of m around 3 
(we will see below it looks like the case) our main concern about the fundamental mode results summarized in Fig. || 
refers to the feasibility of the chosen escape efficiencies. The ratio 7 C /(7 C + 7s) is difficult to maximize in a resonant 
mode because, by its own resonant nature, 7 C must be rather low. Thus, in |2j| it was only of 0.52, while in [Q it 
was 0.36. Even in fl, a doubly resonant system specifically designed to squeeze the fundamental mode, the escape 
efficiency was around 0.9, limiting the maximum squeezing achievable to 90% (in practice, a 52% of noise reduction 
was reached). It appears, then, that nowadays the 77 = 0.99 should be taken rather as an ideal illustrative case. 

In contrast, the ultimate limit for the noise suppression in the harmonic mode (Eq. (5.25)) is pushed up by 
the fundamental mode photon number, opening a way to bypass the usual untouchable limit imposed by the escape 
efficiency of the cavity (as in the fundamental mode). Therefore, the squeezing in the harmonic mode can be arbitrarily 
large under the ideal assumption that the energy load inside the cavity can be also arbitrarily large. However, this is 
not totally true as the model does not take into account the losses in the harmonic mode which necessarily limit the 
degree of noise suppression. We can estimate this limitation assuming the absorption in one single pass through the 
nonlinear material equivalent to the effect of a beam splitter with the adequate reflectivity. Taking an absorption of 
0.6%/cm as in || and a length of 1 cm, the equivalent reflectivity would be of 6 10~ 3 . The spectrum after the beam 
splitter is given by S ou t = 1 + T : Si„ :. Setting : Sj n := — 1 and T = 1 — i?, the ultimate squeezing achievable is 
precisely R = 6 10 -3 , i.e. - 22 dB. In other words, the chosen value of m = 50 in Figs. ^| and || represents more or 
less the maximum the model can stand without the inclusion of the harmonic mode losses. 

Of course, we still cannot assume m — 50 as a realistic limit for the state of the art devices as m depends not 
only on the intra-cavity photon number but on the ratio 1^/7 between the nonlinearity and losses. This ratio must 
be high enough in order to prevent a degradation of the nonlinear optical response in the system as commented 
in the introduction. Besides, this ratio scales down the power available in the external sources. In view of these 
complications, probably the most reliable way of setting the physical scale of m is to compare the results w ith the 
reported experiments. In Q the quoted noise reduction was of -5.2 dB. Setting to zero Afc and /3 in in Eq. ( |5.19b| ) 



15 




correspon ding t o phase-matched SHG, a -5.2 dB squeezing results at m = 2.5, far from the m = 50 limit. Fortunately, 
the limit (5.25) grows up quite quickly for low m's (see Fig. |). Thus, 10 dB of noise suppression are reached at 



m = 5, no such unthinkable value. However, -15 dB of noise reduction requires m — 15, while a -20 dB figure is at the 
m = 50 limit, an order of magnitude higher. New nonlinear materials seem the only possibility for such high squeezing 
degrees. A promising via consists in the use of resonant nonlinearities in asymmetric quantum wells (AQW). Huge 
nonlinearities have been demonstrated in frequency doubling experiments and even a tuning of the nonlinearity with 
a d.c. field j34j|. Obviously, also the absorption is enhanced by the resonance. This can be a problem as the ratio 
v/'y could at the end of the day not be increased. To asses this possibility requires quite a detailed analysis out of 
the scope of the present work. We can foreseen, however, a promising advantage in the fact that the losses in the 
harmonic mode have little influence on the performance. By maintaining a strong two photon resonance but relaxing 
the one photon counterpart (tuning with a d.c. field or by an adequate energy level engineering), the nonlinearity 
would be certainly enhanced while the losses at the fundamental mode would not increase so strongly, thus enhancing 
v/'y. With only one passage through the cavity of the harmonic mode and taking into account that a very thin layer 
of material is capable of SHG [34|] , the corresponding deleterious effect cannot be very large. Even a more exciting 
possibility comes from the recent experimental demonstrations of absorption inhibition in AQW induced by quantum 
interference [|35| — |3T|] . The absorption transparency and the resonant enhancement can be combined using an adequate 
quantum well engineering leading to very efficient frequency doublers (see psfl , where precisely a scheme only resonant 
at the harmonic mode is proposed). 

These are certainly promising perspectives but we should not dismiss the improvements arising at the range of the 
present nonlinear crystals performances. Let us center then, around m = 2.5. As shown in the previous section the 
best strategy corresponds to drive both modes with relative phases ip — <p/2 = tt/2 (negative r]j n ) a nd A = 0. 



In Fig. [10|the noise behavior till m = 5 is displayed for various "distances" to the instability (5.28). Even at half the 
instability r/i n value, the squeezing at m = 2.5 grows from -5.1 dB (69%) to -7.2 dB (80%). The excess noise, on the 
other hand, rapidly increases at low m's but it also saturates quickly to bearable values. The improvement, although 
nothing spectacular is quite substantial. In p7fl it was not reckoned so because the noise suppression was studied as 
function of the input power. Given its nonlinear relation with m the improvement is much slower with respect to this 



variable. Besides the squeezing, the output power is also enhanced. Taking a negative r\i n in Eq. (5.29), the output 
power results in 

Pout oc (2m + |?7m|) 2 , (6.1) 

and thus, the harmonic mode input contributes constructively to it. As shown in Fig. [ll] at half a way of the 
instability the power is nearly doubled. Although from the theoretical point of view the injection of a coherent signal 
in the harmonic mode looks quite harmless, the experimental implementation is not trivial. However, the remarkable 
achievements in Jl0],[ll) with the OPA strongly support the feasibility of the idea. 

Finally, a word of caution about the design of the device. It is important to avoid the setting of oscillations out of 
the fundamental mode (the so called subharmonic pumped OPO p9| , ^0| ), something capable of destroying the noise 
reduction UW. At a first glance, finite values of r]i n would favor the effect by promoting the down conversion. But it 
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FIG. 11. The harmonic output power corresponding to the cases of Fig. 



is not necessarily so as the down conversion is encouraged only for a given range of the relative phase between the 
two driving fields. Thus, for the negative r\i n case studied above, being the harmonic output power maximized, the 
down conversion is minimized. 

To conclude, let us summarize the most relevant results. Firstly, for any system with only one effective mode we 
have given a systematic approach capable of isolating the processing of quantum noise by the dynamic response of the 
system. This dynamic processing is maximized at zero frequency, zero generalized nonlinear detuning (A as defined in 
section [v]) and at a dynamic instability. The static contributions to the noise coming from the different noisy inputs 
can in some cases, move the overall optimum working point away from that corresponding to maximum dynamic 
noise suppression. In spite of such, to have a rule to maximize the dynamic quantum noise suppression resulted very 
useful to characterize the squeezing behavior when applied to a specific optical system. In particular, for the case 
of a singly resonant second order nonlinear device, the squeezing at the fundamental mode is limited by the escape 
efficiency of the cavity, the best working point being within figures of merit of conventional nonlinear crystals. In the 
harmonic mode high squeezing requires new materials but it is only limited by the losses in the non-resonant harmonic 
mode opening the possibility of using multiple quantum wells with resonantly enhanced nonlincaritics. However, with 
standard nonlinear crystals still is possible a substantial improvement with respect to the reported experiments by 
injecting a coherent driving in the harmonic mode. Besides, the output power is highly enhanced. 
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